How much entropy is produced in strongly coupled 
Quark-Gluon Plasma (sQGP) by dissipative effects? 
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We argue that estimates of dissipative effects based on the first-order hydrodynamics with shear 
viscosity are potentially misleading because higher order terms in the gradient expansion of the 
dissipative part of the stress tensor tend to reduce them. Using recently obtained sound dispersion 
relation in thermal 7V=4 supersymmetric plasma, we calculate the resummed effect of these high 
order terms for Bjorken expansion appropriate to RH1C/LHC collisions. A reduction of entropy 
production is found to be substantial, up to an order of magnitude. 

PACS numbers: 



Hydrodynamical description of matter created in high 
energy collisions have been proposed by Landau [l[ more 
than 50 years ago, motivated by large coupling at small 
distance, as followed from the beta functions of QED and 
scalar theories known at the time. Hadronic matter is of 
course described by QCD, in which the coupling runs in 
the opposite way. And yet, recent RHIC experiments 
have shown spectacular collective flows, well described 
by relativistic hydrodynamics. More specifically, one ob- 
served three types of flow: (i) outward expansion in trans- 
verse plane, or radial flow, (ii) azimuthal asymmetry or 
"elliptic flow" 0, [H, as well as recently proposed (iii) 
"conical flow" from quenched jets Q. These observation 
lead to conclusion that QGP at RHfC is a near-perfect 
liquid, in a strongly coupled regime [5]. The issue we 
discuss below is at what "initial time" To one is able to 
start hydrodynamical description of heavy ion collisions, 
without phenomenological/theoretical contradictions. 

Phenomenologically, it was argued in 0, [1[ that elliptic 
flow is especially sensitive to To. Indeed, ballistic motion 
of partons may quickly erase the initial spatial anisotropy 
on which this effect is based. In practice, hydrodynamics 
at RHIC is usually used starting from time To ~ 1/2 fm, 
otherwise the observed ellipticity is not reproduced. 

Can one actually use hydrodynamics reliably at such 
short time? How large is To compared to a relevant "mi- 
croscopic scales" of sQGP? How much dissipation occurs 
in the system at this time? As a measure of that, we will 
calculate below the ratio of the amount of entropy pro- 
duced at t > t to its "primordial" value at t , AS/ So- 

To set up the problem, let us start with a very crude 
dimensional estimate. If we think that the QCD effective 
coupling is large a s ~ 1 and the only reasonable micro- 
scopic length is given by temperature [14j ] , then the rele- 
vant micro-to-macro ratio of scales is simply To To. With 
T ~ 400 Me V at RHIC, one finds this ratio to be close 
to one. We are then lead to a pessimistic conclusion: at 
such time application of any macroscopic theory, thermo- 
or hydro-dynamics, seems to be impossible, since order 
one corrections are expected. 

Let us then do the first approximation, including the 
explicit viscosity term to the first order. Zeroth order (in 



mean free path) stress tensor used in the ideal hydrody- 
namics has the form 



T aJ = ( e + P) u » u v + P9fiv 



(1) 



while dissipative corrections are induced by gradients of 
the velocity field. The well known first order corrections 
are due to shear (77) and bulk (£) viscosities 

ffTW - riiV^Uv + V„u M - ^A^X7 p u p ) + £,(A^X7 pUp ){2) 

In this equation the following projection operator onto 
the matter rest frame was used: 



(3) 



The energy-momentum conservation T pv at this order 
corresponds to Navier-Stokes equation. 

Because colliding nuclei are Lorentz-compressed, the 
largest gradients at early time are longitudinal, along the 
beam direction. The expansion at this time can be ap- 
proximated by well known Bjorken rapidity-independent 
setup Q , in which hydrodynamical equations depend on 
only one coordinate - proper time r — \Jt 2 — x 1 . 
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where we have introduced the entropy density s = (e + 
p)/T. Note that for traceless T plJ (conformally invariant 
plasma), the bulk viscosity £ = 0. 

For reasons which will become clear soon, let us com- 
pare this eqn to another problem, in which large longi- 
tudinal gradients appear as well, namely sound wave in 
the medium. The dispersion relation (the pole position) 
for a sound wave with frequency u> and wave vector q is, 
at small q 
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Notice that the right hand side of (0]) contains precisely 
the same combination of viscosity and thermodynamical 
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parameters as appears in the sound attenuation problem: 
the length T s , which measures directly the magnitude 
of the dissipative corrections. At proper times r ~ T s 
one has to abandon the hydrodynamics altogether, as 
the dissipative corrections cannot be ignored. 

For the entropy production (j4|) the first correction to 
the ideal case is (f — T s /t). Since the correction to one is 
negative, it reduces the rate of the entropy decrease with 
time. Equivalently statement is that the total positive 
sign shows that some amount of entropy is generated 
by the dissipative term. Danielewicz and Gyulassy Q 
have analyzed eq. ^ in great details considering vari- 
ous values of rj. Their results indicate that the entropy 
production can be substantial. 

Our present study is motivated by the following ar- 
gument. If the hydrodynamical description is forced to 
begin at early time tq which is not large compared to 
the intrinsic micro scale 1/T, then limiting dissipative 

effects to the first gradient only (ST^J ) is parametrically 
not justified and higher order terms have to be accounted 
for. Ideally those effects need to be resummed. As a first 
step, however, we may attempt to guess their sign and 
estimate the magnitude. 

Formally one can think of the dissipative part of the 
stress tensor 8T^ U as expended in a series containing all 
derivatives of the velocity field u, 5T^ V being the first 
term in the expansion. In general 3+1 dimensional case 
there are many structures, each entering with a new and 
independent viscosity coefficient. We call them "higher 
order viscosities" and the expansion is somewhat similar 
to twist expansion. For 1+1 Bjorken problem, the ap- 
pearance of the extra terms modifies eq. (U]) , which can 
be written as a series in inverse proper time 
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We have put T here simply for dimensional reasons: 
clearly TV is a micro-to-macro scale ratio which deter- 
mines convergence of these series and the total amount of 
produced entropy. Similarly, the sound wave dispersion 
relation becomes nonlinear as we go beyond the lowest 
order: 



w = K[w(?)] + i$t[w(q)] 
3t[w] q 



(7) 
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Based on T-parity arguments we keep only odd (even) 
powers of q for the real (imaginary) parts of uj. The co- 
efficients c„, r n and rj n are related since they originate 
from the very same gradient expansion of . Although 
both the entropy production series above and sound ab- 
sorption should converge to sign-definite answer, the co- 
efficients of the series may well be of alternating sign (as 
we will see shortly). 



Clearly, keeping these next order terms can be use- 
ful only provided there is some microscopic theory which 
would make it possible to determine the values of the 
high order viscosities. For strongly coupled QCD plasma 
this information is at the moment beyond current theo- 
retical reach, and we have to rely on models. A particu- 
larly useful and widely studied model of QCD plasma is 
Af = 4 supersymmetric plasma, which is also conformal 
(CFT). The AdS/CFT correspondence @] (see @ for re- 
view) relates the strongly coupled gauge theory descrip- 
tion to weakly coupled gravity problem in the background 
of AdSs black hole metric. Remarkably, certain informa- 
tion on higher order viscosities in the CFT plasma can be 
read of from the literature and we exploit this possibility 
below. 

The viscosity-to-entropy ratio (rj/s = 1/4tt) deduced 
from AdS [l0( turns out to be quite a reasonable ap- 
proximation to the values appropriate for the RHIC data 
description. Thus one may hope that the information on 
the higher viscosities gained from the very same model 
can be well trusted as a model for QCD. Admittedly hav- 
ing no convincing argument in favor, we simply assume 
that the viscosity expansion of the QCD plasma displays 
very similar behavior, both qualitative and quantitative, 
as its CFT sister. 

Our estimates are based on the analysis of the quasi- 
normal modes in the AdS black hole background due to 
Kovtun and Starinets [ll|. The dispersion relation for 
the sound mode, calculated in ref. [11(, is shown in FigfT] 
The real and imaginary parts of u> correspond to the ex- 
pressions given in Q. At q — > they agree with the 
leading order hydrodynamical dispersion relation ([5]). 

The first important observation is that the next order 
coefficient 772 is negative, reducing the effect of the first 
one when gradients are large. The second is that 
has maximum at q/2nT ~ 1, and at large q the imaginary 
part starts to decrease. This means that the expansion 
([7]) has a radius of convergence q/2irT ~ 1. 

In order to estimate the effect of higher viscosities on 
the entropy production in the Bjorken setup we first iden- 
tify r in © with 27r/g in ([7]). Second we identify the coef- 
ficients c n with rj n . Both sound attenuation and entropy 
production in question are one dimensional problems as- 
sociated with the same longitudinal gradients and pre- 
sumably the same physics. In practice we use the curve 
for the imaginary part of u> (Fig. [I} as an input for the 
right hand side of ([6|). 

The numerical results are shown in Figs. [5] and [3] in 
which we compare our estimates with the "conventional" 
shear viscosity results from (0}. To be fully consistent 
with the model we set n/s = We also set the initial 

temperature To = 300 MeV while the standard equation 
of state s = 4 ksB T 3 . For the coefficient k$B we use the 
"QCD" value 

kss = ^ Un c - l) 2 + \N c n f \ ■ n } =3; N c = 3 
Fig. [2] presents the results for entropy production as a 



FIG. 2: Entropy production as a function of proper time for 
initial time to = 0.2 fm (left) and to = 0.5 fm (right). The 
initial temperature To = 300 MeV. The dashed (blue) curves 
correspond to the first order (shear) viscosity approximation 
Eq.(U]). The solid curve (red) is the all order dissipative re- 
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FIG. 1: Sound dispersion (real and imaginary parts) obtained 
from the analysis of quasinormal modes in the AdS black hole 
background. The result and figure are taken from Ref.jlll]. 



FIG. 3: Fraction of entropy produced during the hydro phase 
as a function of initial proper time. The initial temperature 
To = 300 MeV. The left (blue) points correspond to the first 
order (shear) viscosity approximation. The right (red) points 
are for the all order resummation. 



function of proper time for two initial times t = 0.2 fm 
and r = 0.5 fm. The dashed lines correspond to the first 
order result ^ while the solid curves include the higher 
order viscosity corrections. Noticeably there is a dra- 
matic effect toward reduction of the entropy production 
as we start the hydro evolution at earlier times (the ef- 
fect is almost invisible on the temperature profile). This 
is the central message of the present paper. 

Fig. [3] illustrates the relative amount of entropy pro- 
duced during the hydro phase as a function of initial time. 
If the fist order hydrodynamics is launched at very early 
times, the hydro phase produces too large amount of en- 
tropy, up to 250%. (Such a large discrepancy is not seen 
in the RHIC data.) In sharp contrast, the results from 
the resummed viscous hydrodynamics is very stable, and 
does not produce more than some 25% of initial entropy, 
even if pushed to start from extremely early times. The 
right figure displays the absence of any pathological ex- 
plosion at small to. 

It is worth commenting that we carried the analysis 
using the minimal value for the ratio r]/s = 1/47T. We 
expect that if this ratio is taken larger, the discrepancy 
between the first order dissipative hydro and all orders 
will be even stronger. 

Before concluding this paper we note that a practi- 
cal implementation of relativistic viscous hydrodynamics 
had followed Israel-Stewart second order formalism (for 
recent publications see [l2l |) in which one introduces ad- 



ditional parameter , the relaxation time for the system. 
Then the dissipative part of the stress tensor is found as 
a solution of an evolution equation, with the relaxation 
time being its parameter. For the Bjorken setup, the dis- 
sipative tensor thus obtained has all powers in 1/t and 
might resemble the expansion in ^ and ([7]). The use of 
AdS/CFT may shed light on the interrelation between 
the two approaches: the first step in this direction has 
been made recently [l3j |. resulting in numerically very 
small relaxation time. 

Finally, why can it be that macroscopic approaches like 
hydrodynamics can be rather accurate at such a short 
time scale? Trying to answer this central question one 
should keep in mind that 1/T is not the shortest micro- 
scopic scale. The inter-parton distance is much smaller, 

1 /3 

~ 1/(T * N do f) where the number of effective degrees of 
freedom N do f — 40 in QCD while N dof — N% — > oo in 
the AdS/CFT approach. 

In summary, we have argued that the higher order dis- 
sipative terms strongly reduce the effect of the usual vis- 
cosity. Therefore an "effective" viscosity-to-entropy ratio 
found from comparison Navier-Stokes results to experi- 
ment, can even be below the (proposed) lower bound of 
1/47T. We conclude that it is not impossible to use a hy- 
drodynamic description of RHIC collision starting from 
very early times. In particular, our study suggests that 
the final entropy observed and its "primordial" value ob- 
tained right after collision should indeed match, with an 
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accuracy of 10-20 percent. 
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